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I. INTRODUCTION 

The past decade has witnessed great advances in both 
the fields of cold atoms and cavity quantum electrody- 
namics (Cavity-QED), and the overlap between the two 
fields is ever-growing. A remarkable achievement in this 
direction is the successful coupling of a Bosc-Einstein 
condensate to a quantized field mode of a high-finesse 
optical cavity [J H|. Besides that, deterministic loading 
of individual atoms in a micro-cavity is demonstrated [3| 
and submicronpositioning of single atoms in the cavity 
is achieved [H, |j| , which allows control of the atom- field 
coupling via its position-dependence. 

Theoretically, Mekhov, Maschler, and Ritsch proposed 
to probe the superfluid-insulator transition of cold atoms 
in optical lattices by the transmission spectra of an op- 
tical cavity ||. The atoms couple to a quantized cavity 
mode dispersively and hence act as some moving refrac- 
tive media in the cavity. The cavity transmission spectra 
directly reflects the quantum or classical distributions of 
the atoms, which characterize the superfluid or insula- 
tor phases respectively. This non-destructive proposal 
exploits the fact that in the domain of strong coupling, 
even one atom is enough to shift the cavity resonance sig- 
nificantly. Techniques based on this knowledge have been 
developed to detect the existence of atoms in a cavity Q , 
and most recently, been employed to study the correla- 
tion, statistics and dynamics of matter- wave fields [7[ . 

From the point of view of atomic optics and quantum 
information, Ref. @ also provides us with a model of rich 
coupled atom-field dynamics ||. The atoms effectively 
influence the field dynamics by shifting the resonance of 
the field mode, while in turn the field intensity determines 
the dipole potential for the atoms. The former effect is 
essential for the result of Ref. [f| and is treated in detail. 
However, the latter effect is neglected. The atomic dy- 
namics is avoided by prescribing a state (phase) for the 
atoms. Furthermore, the interaction and the coupling to 
the environment may induce entanglement between the 
atomic and field subsystems, and cause decoherence of 
the subsystems, respectively. All these aspects of the 
system are rarely investigated in Ref. [j| . 



The purpose of this paper is to investigate the dynam- 
ics of the composite atom- field system, with the emphasis 
on the interplay between the two sides, the correlation 
and the entanglement between them. We shall consider 
a "two-site version" of the model presented in Ref. [f|. 
Atoms are trapped in a double-well potential and inter- 
act dispersively with a damped and driven field mode. 
The two traps are placed asymmetric to the field mode 
so that the atomic tunneling dynamics is coupled to the 
field dynamics. Under the two-mode approximation, the 
freedoms of the atoms are reduced to minimum and can 
be taken into full account. To gain insight into the dy- 
namics of the system, we assume the system starts from 
an initial state and evolves towards the steady state. We 
find that this process involves two distinct time scales, 
one is the atomic tunneling rate and the other the cavity 
loss rate, with the latter much faster than the former. 
These two incommensurate time scales lead to distinct 
temporal structures of the dynamics. In the short-time 
interval, where the atomic tunneling can be neglected, it 
is found that the model we consider is analogous to the 
Dicke model in the dispersive regime, for which a good 
understanding exists 9]. Detailed analytical results are 
obtained and, by the way, the main result of Ref. [f| is 
recovered. In the long-time interval, however, the atomic 
tunneling plays an important role. Strong population 
transfer between different atomic states is observed, and 
the system displays substantially different behavior than 
in the short-time interval. 



This paper is organized as follows. In Sec.|TT]the basic 
model is introduced and the Hamiltonian of the atom- 
field system is derived under the two-mode approxima- 
tion. Then in Sec. IIII[ based on the master equation, 
the short-time and long-time behaviors are investigated 
both analytically and numerically. Finally, our results 
are summarized in Sec. IIVI The connection with the 
Dicke model is discussed in Appendix A and some useful 
inequalities are derived in Appendix B. 
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II. BASIC MODEL AND THE HAMILTONIAN 

In this work, we consider the combination of a double- 
well and an optical cavity, two paradigm models in 
physics. We assume N two- level bosonic atoms with mass 
m and transition frequency u) a are trapped in a double- 
well potential V(x) and loaded in an optical cavity, where 
they interact with a cavity field mode with frequency to c . 
The cavity is coherently pumped through the mirror by 
a weak laser with frequency io v and amplitude ry. We 
also assume the atom-field detuning is much larger than 
the atomic spontaneous emission rate and the Rabi fre- 
quency. Under this condition, the atomic upper level can 
be adiabatically eliminated [l(J EH, i- e -> the atomic in- 
ternal dynamics is neglected. 

After adiabatic elimination of the atomic upper state, 
the single-atom-plus-field Hamiltonian in the frame ro- 
tating at the frequency of the pumping field is [Icl |: 



H o — H. 



ph 



(1) 



where H p h is the rotating frame Hamiltonian for the 
driven field. 



H ph = —Aa'a + r](a + a^), 



(2) 



with A = uj p — lu c being the pump-cavity detuning, and 
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H s = ^- + V(x) + u 2 (x)(U Q a*a), (3) 

which is the Hamiltonian for a single atom in the super- 
position of the classical potential V(x) and the quantum 
potential u 2 (x)U a^a Here u(x) is the field mode 

function with its magnitude at the antinode normalized 
to unity. The parameter Uq = g^jiyjc — u a ), with go 
being the atom-field coupling at the antinode. 

The many-atom-plus-field Hamiltonian, taking into ac- 
count the direct interaction between the atoms which is 
characterized by the s-wave scattering length a s , is: 



H = H ph + J d 3 x^(x)H a y(x) 



1 Ana. 



2 m 



<Px&(x)&(x)9(x)V(x), (4) 



where ^(x) is the atomic field operator and we take Ti = 1 
here and henceforth. Under the two-mode approximation 
for the atomic freedoms, the validity of which has been 
well established Ref. [l3| , the atomic field operator 
ty(x) has two contributions: 



^(x) = biwi(x) + b 2 w 2 (x) 

= b\w(x — x\) + b 2 w(x 



X 2 



(5) 



Here we assume a symmetric double well with the two 
minima at x\ and x 2 ■ The two modes ^1(2;) and w 2 (x) 
are localized in the left and right traps respectively, and 
satisfy the orthonormal relation J d 3 xw* (x)uij(x) = Sij, 



= 1, 2. The operator b\ (pi) (i — 1, 2) creates (anni- 
hilates) an atom in the mode Wi (x) . Substituting Eq. ([5]) 
into Eq. ([4]), and keeping only terms with dominating 
contributions, we obtain 

H = H p h + H a + Hi n t, (6) 

where H a is the Hamiltonian for the atomic subsystem, 
H a = -t(b\b 2 + b\bx) + ^{nx{n x - 1) + n 2 {n 2 - 1)). (7) 

Here we introduce the atom number operators rii = b-bi 
(i = 1,2) and drop the term associated with the zero- 
point energy. The atomic tunneling rate t and the on-site 
interaction energy u are defined as 

-t = J d z xwl{x)(-^- + V{x))w 2 {x), (8) 
47ra„ 



2m 

d 3 x\w h2 (x)\ 4 . 



(9) 



Hi n t describes the interaction between the atoms and the 
field, 

H lnt = J d 3 x^(x)u 2 (x)^(x)(U a^a) 

~ (Jim + J 2 n 2 )(U Q a^a), (10) 
where the dimensionless coefficients J\ 2 are defined as 



■h. 



d 3 xu 2 (x)\w 1<2 (x)\'' 



(11) 



which reflect the overlap between the atomic modes and 
the field mode. Note that J\ j2 are bounded, < J\_ 2 < 1, 
which follows from the normalization conditions of u(x) 
and ^1,2(3;). If the field mode u(x) varies slowly in the 
range of the spread of the atomic modeSj we can take the 
"tight confinement approximation" 0, [lij Ji ~ u 2 (xi) 
(i = 1,2). It is clear from Eq. (fTQ)) that the interac- 
tion between the atoms and the field is twofold. For the 
atoms, the depths of the two traps are shifted while for 
the field the energy per photon is renormalized. 

We shall discriminate two different cases: (i) J\ = J 2 \ 
(ii) J\ ^ J 2 . The former case is trivial, because in this 
case the dynamics of the atoms and the field is decoupled, 
the field is indifferent to the distribution of the atoms in 
the two traps. Thus we concentrate on the case J\ 7^ J 2 . 
Without loss of generality, we assume J\ = 1, J2 = 0. 
This is always reasonable mathematically because we can 
define two effective parameters, A' = A — UqJ 2 N, Uq = 
Uq(Ji — J2), and rewrite the Hamiltonian as 

H = -t(b\b 2 + blh) + |(ni(ni - 1) + n 2 (n 2 - 1)) 



+UQa)ani — A'a^a + rj(a + a 1 ), 



(12) 



then effectively we have J\ = 1, J2 = 0. Experimentally, 
excellent control of the position of a single atom relative 
to the cavity mode has been demonstrated, so atom-field 
coupling can be tailored as wanted [1, 0, [l5|, [l(| • In the 
following, we shall omit the prime for notational simplic- 
ity. 
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III. ANALYTICAL AND NUMERICAL 
ANALYSIS BASED ON MASTER EQUATION 

The Hamiltonian derived above controls the coherent 
evolution of the atom-field system. However, we still have 
to take the dissipation into account, which comes from 
the cavity loss in the model we consider. The overall evo- 
lution of the system is governed by the master equation: 

p = —i[H, p] + n{2apa) — a)ap — pa) a) = Cp. (13) 

Here p is the density matrix of the atom-field system in 
the rotating frame, and n is the cavity loss rate. Note 
that generally the frequency of the cavity mode falls in 
the optical regime, hence the environment can be treated 
as at zero temperature. The master equation will be our 
starting point for the rest of the paper. 

As for the dynamics of our system, we stress that there 
are two distinct time scales . One is the inverse of the 
atomic tunneling rate t _1 , the other being that of the 
cavity loss rate ft -1 . They are the characteristic times of 
the atomic and field subsystems, respectively. In typical 
experimental situations, k is of order 10 6 Hz, while t (and 
u) is of order 10 3 Hz at most [13j. This means that gener- 
ally there is a hierarchy i -1 ^> kT 1 . The identification of 
two different time scales leads us to classify the dynamics 
of the system into short- and long- time behaviors, which 
correspond to two disjoint time intervals, (i) < r <C t 
and (ii) t 3> t , respectively. In the short-time inter- 
val, the atomic tunneling is "frozen" . However, we still 
expect the system to display some non-trivial behaviors, 
because this time may be long in unit of In the 

long-time interval, the atomic tunneling may eventually 
give rise to some important results and should be taken 
into full account. Specifically we divide the Hamiltonian 
H into tunneling and non-tunneling terms, 

H = H t + H non , (14) 

with 

H t = -t(blb 2 + blh), (15) 

H non = —AaJa + r)(a + a') + Uoa'ani 
u 

+ ^{ni(n 1 -l) + n 2 (n 2 -l)), (16) 

and rewrite the master equation as 

p= -i[H non , p] + n(2apa) -a f ap- pa) a)-i[H t , p\. (17) 

The last term will be neglected (kept) in the short- (long- 
) time intervals, respectively. In the following we shall 
investigate the behavior of the system in the two time 
intervals both analytically and numerically. 

A. Short-time behavior 

Let us assume initially the atoms are in the ground 
state \G) of the Hamiltonian H a , while the field is in the 



vacuum state |0)/, i.e., 

p(0) = \G)(G\®\0) f f(Q\. (18) 

Then at r = the pump is turned on and the system 
evolves according to the master equation (fT3|) . In gen- 
eral, solving a master equation analytically exactly is a 
formidable task, so we will resort to numerical methods 
as we do. However, in the short-time interval, as men- 
tioned above, we may neglect the tunneling term and 
approximate the master equation by 

p = -i[H non , p] + K,(2apa) - a) ap— pa) a) = C non p. (19) 

As pointed out in Appendix A, H non can be mapped into 
the Dicke model in the dispersive regime, up to some 
minor differences. The dynamics of the Dicke model in a 
driven and damped cavity, in the dispersive regime, has 
been studied in detail in Rcf. [9]. Here we shall follow 
the techniques there. 

Under the transformation to another reference frame, 

~ iAr —iAr /on\ 

p = e pe , (20) 

u 

A = -(n 1 (n 1 -l)+n 2 (n 2 -l)), (21) 
the master equation fTS")) takes the form 

p = —i{H, p] + n(2apa) — a) ap — pa) a), (22) 
with the simplified Hamiltonian 

H = -Aa f a + rj(a + a)) + U a)ani. (23) 

Note that H is diagonal in the atomic space. This leads 
us to expand the density matrix p as 

N 

p= ^2 \m){n\<8)p mn , (24) 

m,n— 

where |m) = \m,N — m) denotes the atomic state with 
m atoms in the left trap and (N — m) atoms in the ritht 
trap, and p mn = (m\p\n), which is still an operater in 
the field space. The initial values of the p m „'s are 

p mn (0) = (m\G)(G\n) ■ \0) ff (0\. (25) 

In terms of p mn , the atomic and field density operators 
are respectively, 

N 

Pa = trf(p) = ^ tr f (p mn )\m)(n\, (26) 

m,n— 
N 

Pf = tr a (p) = ^ Pmm- (27) 
m=0 

It is straightforward to obtain the time evolution equa- 
tion of the operators p mn from the master equation (|22[) , 

p mn = -i[-(A-U p)a)a + r)(a + a)),p mn ] 
+n(2ap mn a) - a) ap mn - p mn a) a) 
-iU q( & O'Pmn ~h Pmn a) a), (28) 
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where p = (m + n)/2, q = (to — n)/2. The general 
solution of this equation is derived in Ref. [§] by applying 
the dynamical symmetry method. The result is thorough 
but complicated, so we will cite it only when we have to. 

For the diagonal cases with to = n, q = 0, Eq. ([28]) 
reduces to the master equation describing the dynamics 
of a single field mode subjected to damping and pumping. 
Up to a constant coefficient, we make the ansatz that the 
solution of Eq. ([28]) in this case takes the form 



Pmm (^~) 



l ( T )>//("m(T)|, 



(29) 



where |a m (r))/ denotes the field coherent state, and 
a m (0) = 0. Substituting ([29]) into Eq. ([28]). after some 
operator manipulations, we find that the equation is sat- 
isfied if 



a m (r) = -(« - i(A - U m))a m (T) - it]. 
This equation is readily solved by 

a m (r) = a m (oo)(l - e -(^(A-c/ m))r ); 

with 

— IT] 



a m (oo) = 



k — i(A — U$m) 



(30) 
(31) 
(32) 



Recalling Eqs. (|25|) and (27), we have the field density 
operator 



N 



PS 



(r)=J2\HG)\ 2 \^m(r)) ff (a m (r)\, (33) 



m=0 



which is an incoherent superposition of a series of co- 
herent states. Due to the atom-field coupling, both the 
weights of the coherent states and the coherent states 
themselves, depend on the initial atomic state. For times 
r ^> a m (t) saturates to the value a m (oo) and the 
field approaches the quasi-steady state 



Pf(°°) 



N 

E 



(m|G)| 2 |a m (oo)) // (a m (oo)|. (34) 



The average photon number in the quasi-steady state is 



N 



Ja)=J2\(m\G)f 



m=0 



k 2 + (A - U my 



(35) 



Note that the field approaches its quasi-steady state in 
a time of order which is well within the short-time 
interval < r < f 1 . This indicates that the analysis 
above is self-consistent. We refer to this "steady state" 
of the field as quasi-steady state so as to differentiate it 
from the true steady state in the long-time interval. 

For off-diagonal cases with to ^ n, the last term in 
Eq. ([2"8"]) is nonzero. As pointed out in Ref. this 
non-unitary term will result in the complete disappear- 
ance of the operators p mn , that is, the complete coher- 
ence loss of the atomic subsystem. Explicitly, 



mn 
\Pa 



\trf(p mn )\ oc exp(-T/T mn ), 



(36) 



A 

+ ro 

+ ro 
v 




FIG. 1: (color online). Normalized photon number 
(a) a) / (rf / k 2 ) as a function of the pump-cavity detuning A, 
with the master equation cut off at four different times. The 
quasi-steady state result (solid line) is shown for comparison. 
The two lines corresponding to kt = (20, 200) have been up 
shifted 0.15 and 0.30 respectively, unless they coincide with 
the solid line. The parameters are (t, u) = 2tv x (400, 200)Hz, 
(k, Uo, 77) = 2tt x (1.5, 6.0, 0.1) x 10 6 Hz. The number of atoms 
is N = 2. 



with the (to, n)-dependent characteristic time 

= [k 2 + (A - U m) 2 }[n 2 + (A - U n) 2 
KU^r\ 2 (m — n) 2 



(37) 



If all the r mn 's are much smaller than t , then even- 
tually the atomic subsystem will reach a purely mixed 
state, 



N 



p a (oo) = J2 \{m\G)\ 2 \m){m\ 



(38) 



m=0 

and the atom-field system is merely classically correlated, 

JV 



p(oo) = ^ |(HG)| 2 |to)(to| ® |am(oo))//(a m (oo)|. 



m=0 



(39) 



All the results derived above are based on the approxi- 
mation that the atomic tunneling is negligible in the short 
time interval. The quality of this approximation is well 
demonstrated in Figs. [1] and [2] There we show the time 
evolution of the photon number and off-diagonal element 
p® . The results are obtained by numerically integrating 
the master equation (| 13|) with the atomic tunneling being 
taken into account. As shown in Fig.[TJ within a time of 
order K , the photon number builds up and saturates to 
the value given by Eq. ([55]) . Then it holds on to times of 
order 10 3 /k before signatures of deviation from the ap- 
proximation arise. The excellent agreement between the 
analytical and numerical results is again demonstrated in 
the decay of the off-diagonal element p^ 1 in Fig. [21 



5 



CD 



0.35 



0.30 



0.20 



0.15 



- kt=20, numerical 

■ kt=20, approx 
kt=200, numerical 

■ kt=200, approx 



-4 -3 -2 -1 
A/U„ 

FIG. 2: (color online), 
at two time sections, kt = (20, 200). Analytical approximate 
results according to Eqs. (|36|) and (I37p and numerical results 
based on master equation (|13(l are shown for comparison. The 
parameters are the same as in Fig. [T] 
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FIG. 3: (color online). Normalized photon number at steady 
state {a* a) st / {if / k 2 ) as a function of the pump-cavity de- 
tuning A for three different pump strengthes. The dashed 
line corresponding to the quasi-steady state results given by 
Eq. (|35|l is shown for comparison. The parameters are the 
same as in Fig. [Tj 



B. Long-time behavior 

As shown above, in the short-time interval, the atomic 
tunneling can be neglected. However, in the long-time 
interval, where the system enters the steady state p st , 
the atomic tunneling does play an important role. An 
analytical exact solution of p s t is unavailable, so we rely 
on numerical methods 18]. In Fig. [3] we show the the 
normalized photon number in steady state as a function 
of the detuning A with the pump strength varied. The 
difference between the long-time steady state results and 
the short-time quasi-steady state result is apparent. A 
striking feature of the spectrum in steady state is that the 
peaks are almost of equal height and in particular, in the 
weak pump limit (j]/n <C 1), the height converges to some 
value around 1 /3 (take into account the overlap between 
the peaks). In contrast, for the specific set of parameters 
in our numerical calculations, (t, u) — 2n x (400, 200) 
Hz and N = 2, the quasi-steady state result Eq. (155)) 
predicts the heights of the three peaks to be 0.23, 0.53 
and 0.23, respectively. The difference between the steady 
state and quasi-steady state may be more directly re- 
vealed in Fig. |4l where we present the diagonal element 
and off-diagonal element p^ 1 of the atomic density 
matrix p a as the detuning and pump strength are var- 
ied. From Fig. Q|a) we see that when the detuning is far 
from all possible resonances, the element is around 
1/3 regardless of the pump strength; and in the limit of 
weak pump, is around 1/3 in the whole range of the 
detuning. From Fig.fJJb) we see the off-diagonal element 
p^ 1 is far less than unity in the domain of A and r\ we 
consider. Other diagonal and off-diagonal elements have 
similar behavior and hence are not shown. 

We also investigated the cases with N ^ 2, and some 
common features are found. That is, as long as the con- 
dition (t, u) <c (/c, Uq) is fulfilled, in the weak pump 
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FIG. 4: (color online), (a) Diagonal element p°° and (b) off- 
diagonal element pi 1 of the reduced atomic density matrix p a 
at steady state. The parameters are the same as in Fig. [T] 



limit the normalized photon number in steady state 
(a* a) s t / [r] 2 / k 2 ) as a function of the detuning A is the 
superposition of N + 1 lorentzians, which are centered at 
A = U s (s = 0, 1, . . . N), and of heights nearly 1/(N+1). 
Besides that, the diagonal elements of the atomic density 
matrix converge to values around 1/{N +1), while all 
the off-diagonal elements are vanishingly small, i.e., the 
atomic subsystem is in a nearly absolutely "unpolarized" 
mixed state. 

The analysis of the short-time behaviors can in fact 
help us to understand the features of the steady state. 
The steady state p s t satisfies 



= Hp st = CnonPst - i[H t ,p st ]. 



(40) 



Because t (k, Uo), we shall treat the second term as 
a perturbation over the first term, for which we have 
analytical results. Assume that 



Pst 



p° st 



Pit: 



(41) 
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where p° st is of zeroth order in t/n, while components of 
higher orders in t/ k are included in p st , p° st satisfies the 
equation C non p° st — 0. According to the analysis in the 
previous subsection, its general solution is 

N 

p° st = C m \m)(m\ <g |a m (oo))(a m (oo)|, (42) 

m=0 

with the coefficients C m being arbitrary. Note that p° st 
is diagonal in the atomic space, which implies that the 
off-diagonal elements of the atomic density matrix must 
come from p\ t and hence are at least of order t/ k. This 
explains why the off-diagonal elements are vanishingly 
small as revealed by the numerical calculations. The 
physical picture is that, via the atom- field coupling and 
the dissipation, the coherence of the atomic subsystem 
is greatly depleted, the remaining weak coherence is just 
due to the finite atomic tunneling. 

The knowledge of the off-diagonal elements of the 
atomic density matrix allows us to understand the be- 
havior of the diagonal elements and the photon number, 
at least in the weak pump limit. In steady state, we have 
the following equation for an arbitrary operator O, 

0=(6) st = -t([d,H}) st + K ([a\d}a-a^a,d}) st . (43) 

Let O = \m)(m + l|,0<m<JV — 1, then we obtain 

= /(m + l)((|m)(m|) st -(|m + l)(m + l|) st ) 

+/(m + 2)(\m)(m + 2\) st - /(m)(|m - l)(m+ l\) st 
u 

— -(2m+l-N)(\m)(m + l\) st 

-^(\m)(m + l\Ja) st , (44) 

where f(m) — y/ (m + l)(iV — m). In the weak pump 
limit, the second and third terms on the right hand side 
is of order t/n, the fourth term u/k, while the fifth term 
(tj/k) 2 , so to zeroth order in tj k, u/k and ij/k, we have 

(\m)(m\) st - (\m + l)(m + l|) st = 0. (45) 

This equation, together with the normalization condition 
t r (Pa) = L means that to zeroth order in t/n, u/k and 

7]/K, 

(\m)(m\) st = j^j. (46) 

Returning to Eqs. (|4"Tj) and (|4"2"|) , we see that in the weak 
pump limit (the condition t, u <C k is spontaneously sat- 
isfied), the steady state is well approximated by 

f N 

Pst ~ ivTT ^ l TO ^ m ' ® \ a m(oo)) ff (a m (oo)\. (47) 

m=0 

The photon number (a^a) st in this limit is given by 

I N 2 

(atfl> - = FTT ^ k 2 + (A - U m) 2 ' (48) 

m— ' 



Eqs. (|46| and (|48[) account for the weak pump steady 
state features. In fact the photon number is always di- 
rectly determined the diagonal elements (\m) (m\) s t , not 
limited to the weak pump limit. Let O — a 1 a and 
a\m)(m\ in Eq. (|43[) . we have 

= irj{{a) st - {a)) st ) - 2/«(a t a) st 

N 

= ir lY (( a \ m ){ m \)st - c.c.) -2K{a ] a) st , (49) 

m— 

= it([a\m){m\, b\b 2 + b\b{\) st - irj (\m)(m\) st 

-(K-i(A-U m)){a\m){m\) at , (50) 

where c.c. stands for complex conjugate. As shown in Ap- 
pendix B, the first term on the right hand side of Eq. (50) 
can be safely neglected compared to the second term. By 
neglecting it we get a set of algebraic equations of (a'a) s t, 
(a\m)(m\) st . We solve 

N 2 

(JaU - E (HH)st ■ K 2 + {/ ?_ Uom) 2 > (51) 

jn—0 

which is valid for arbitrary values of r\. 

An important question of concern is what is the time 
scale for the system to approach the steady state. The 
general solution of a master equation like Eq. ([13]) with a 
time-independent Liouvillian can be written as a sum of 
a series of complex exponentials, 

p{r) =^2a j exp{s j T), (52) 

3 

where Sj = —Rj+ilj, (Rj,Ij € TV) are the eigenvalues of 
the Liouvillian, while the coefficients a 3 - are determined 
by the initial conditions. As well known, the Liouvillian 
is singular and has at least one zero eigenvalue which 
correspond(s) to the steady state(s), and all the non- 
zero eigenvalues have negative real parts. This ensures 
that p(r) converges to the steady state(s) in the limit 
of t — > oo. Obviously, the time scale of this process is 
set by the inverse of the least modulus real part of the 
eigenvalues. We define 

T m ax = max{-t-; Rj 0}. (53) 
j Hj 

In FigO we show KT max as a function of the detuning. 
Note that KT max is of order f0 4 ~ 10 6 in the regime 
|A/Z7o| < 5 and it diverges as |A| — * oo. Compared 
with Fig. [1] this clearly demonstrates that the long-time 
and short-time behaviors lie in two well separated time 
intervals. Moreover, it clarifies a point that, when the 
pump is far from all possible resonances, the influence of 
the pump on the system is negligible, so it would take the 
system a time experimentally unaccessible to reach the 
steady state. Of course, the situation is still good in the 
regime < A/J7o < 2, where the time scale is of order 
10ms. At present, individual atoms have been trapped 
and detected in an optical cavity for time scales exceeding 
15s [3], so we expect that the steady state features may 
also have the possibility to be observed in the future. 
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3x1 6 




FIG. 5: (color online). Timescale T max in unit of 1/k for the 
system approaching the steady state. Inset: close-up of the 
curve in the regime < A/Uo < 2. The parameters are the 
same as in Fig. [1] 



IV. SUMMARY AND REMARKS 

We investigated the dynamics of a dispersively inter- 
acting atom-field system, with the slowly varying atomic 
interwell tunneling coupled with the fast varying field dy- 
namics. Depending on the role of the atomic tunneling, 
the dynamics of the system was classified into short-time 
and long-time behaviors. 

In the short-time interval (0 < r <C t^ 1 ), as verified 
by the numerical calculations, the atomic tunneling can 
be neglected. We recovered the result of Ref. [5| in the 
"two-site" case, and went beyond to obtain a more de- 
tailed picture of the dynamics of the atom-field system, 
such as the decoherence of the atomic subsystem, the 
correlation between the atomic and field subsystems. In 
our analysis, a central observation is the analogy between 
the model we consider and the well known Dicke model 
in the dispersive regime. In fact, many results are di- 
rectly borrowed from previous works on Dicke model Q . 
Of course, we stress that this similarity is not essential. 
It is the dispersive nature of the atom-field coupling that 
counts. As can be seen from our procedures, similar tech- 
niques and results apply also to the many-site cases, e.g., 
the original model in Ref. [f|. 

As for the long-time behavior, we were primarily in- 
terested in the steady state. If the atomic tunneling is 
absent, the steady state of the system is in the form 
of Eq. (|42| . The atomic and field subsystems are only 
classically correlated, and the populations of different 
atomic states are absolutely determined by the initial 
state. However, the presence of atomic tunneling leads 
to strong population transfers between the atomic states. 
A remarkable feature is that, in the weak pump limit, 
the atomic states are almost equally populated, which is 
substantially different from the ground state atomic dis- 
tribution. We also quantitatively investigated the time 
scale of reaching the steady state and found that it lies 
well in the long time interval and is accessible under the 
experimental situations at present. 



Finally, we have some remarks about the experimental 
implementation of the model we discussed in this work. 
The double- well may be constructed by two adjacent op- 
tical dipole traps between which the distance can be ad- 
justed as in Ref. [19j. An important feature of these traps 
is the extremely small focal spot, a beam waist radius of 
wo < 1/im is achieved. This helps to confine the atom in 
a very small volume and validate the "tight confinement 
approximation". The case J\ = 1, J% = occurs when 
one atomic mode is localized in an antinode of the field 
mode, and the other in a node. 
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APPENDIX A: CONNECTION WITH THE 
DICKE MODEL 

In terms of Schwin ger' s representation of the angular 
momentum operators [20| . 



Sy = -(b^-blh), 

S z = \{b\bi-blb 2 ), 
the Hamiltonian H can be rewritten as 



with 



2tS x , 
UqN 
2 

-u(S 2 z 



H = H t + H n 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 



A)a t a + -y ( 2atfl + l)S z + r,(a + a f ) 



AT2 



2 ; 



U 



(A6) 



Up to terms diagonal in the ^-representation, H non cor- 
responds to the Dicke model in the dispersive regime 0] , 
with cavity-pump detuning (^9p- — A), effective atom- 
field coupling ^ and pump strength r\. It is the two 
center-of-mass motion modes that correspond to the two 
atomic internal levels involved in the Dicke model. 

In this formalism, it is clear that the role of Ht is to in- 
duce transitions between different eigenstates of S z (that 
is, the |m)'s, S z \m) = (m — -y)|m)), with amplitudes 
proportional to t. However, since t k, this process 

can be neglected in the short-time interval. 
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APPENDIX B: SOME USEFUL INEQUALITIES 



We first introduce an inequality [21 1 



\(A f B)\ 2 < {A^A)(B^B), 



(Bl) 



where A, B are two arbitrary operators (hermitian or 
non-hermitian) , and the average is taken over an arbi- 
trary state (pure or mixed). Let A = a, B = I (unity 
operator) in Eq. (|Bip , we obtain 



(a f )(a) < (a f a). 



(B2) 



For steady states, we have this inequality plus the con- 
straint Eq. (|49| . hence 

(a f a) st > I[(( a t) st + (a) st ) 2 -((at} st -(a) st ) 2 ] 



i((at) st + <a> st ) 2 + ^(ata) 2 t 



(B3) 



which yields an upper bound for the photon number 

,2 



(a' a) a t < 



T 



(B4) 



This can be understood as the maximum photon number 
occurs when the probe is at resonance with the cavity. 

We then show that it is legitimate to neglect the first 
term on the right hand side of Eq. (f5t!)) . which is 

it([a\m)( m \,b{b 2 + blh}) st 
= itf(m + l)(a\m)(m + l\) st - itf(m + l)(a\m + l)(m\) st 
+itf(m)(a\m){m - l|) st - itf{m){a\m - l)(m|) . (B5) 



It is ready to show that the ratios of the four terms 
to — irj(\m)(m\) s t are at least of order t/n, as long as 
(\m)(m\} s t is of order unity. Let us take the third term 
for an example. Applying inequalities (|B1[) and (|B4[) . 



\itf(m){a\m){m-l\) st \ < tf(m) [(a^a) st {\m){m\) st ] 1/2 



t„\V2 



< tf(m)(a"a 

K 



(B6) 



and similar results hold for the other three terms. This 
guarantees that the first term on the right hand of 
Eq. (f50]) is much less than the second term. 
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